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We report a way of wave function estimation for the density matrix renormalization group 
(DMRG) method apphed to quantum systems, which has 2-site modulation, when the system 
size extension is necessary in both the finite and the infinite algorithms. The estimation is 
performed by renormalization group (RG) transformation applied to the ground state wave 
function, which is represented as the matrix product. This RG scheme is known as the product 
wave function renormalization group (PWFRG) method. In order to treat the 2-site modula- 
tion, the operation of the RG transformation is shifted by amount of 2 lattice sites. It turns out 
that this 2-site shift algorithm provides better wave function estimation in the thermodynamic 
limit, compared with the previously known PWFRG algorithm. 
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1. Introduction 

Variational estimation of minimum eigenvalues of 
quantum Hamiltonians and maximum eigenvalues of 
classical transfer matrices has been investigated as a non 
perturbative way of analysis in condensed matter sys- 
tems. The Kramers- Wannier approximation applied to 
the two-dimensional (2D) Ising model is one of the early 
example.^' ^ Baxter extended this formalism by introduc- 
ing auxiliary variable, and established the way of cor- 
ner transfer matrix. ^^^ In the field of ID quantum spin 
system, the variational estimation of the ground state 
energy by Nightingale and Blote is one of the earliest 
trial. ^ Quantum state constructed as a product of local 
factors occasionally represent exact ground state, or are 
good variational states.^' ^"'^■^^ Such states are known as 
the matrix product state (MPS), or the finitely correlated 
state.^^^^ Practical and flexible use of the MPS for eigen- 
value problems began with the density matrix renormal- 
ization group (DMRG) method, ^"^ which has been ap- 
plied to various problems in low dimensional correlated 
systems.^"*' ^^ The variational structure in DMRG formal- 
ism mediaged by MPS is revealed by Ostlund and Rom- 
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It is known that numerical calculation in DMRG 
method can be accelerated by explicit use of the ma- 
trix product structure of the variational state, especially 
when the method is applied to finite size ID quan- 
tum systems. ^^'^^ This acceleration procedure can be 
regarded as renormalization group (RG) transformation 
applied to the ground state wave function. It is possi- 
ble to introduce this way of RG transformation to the 
infinite system DMRG method, where the acceleration 
procedure is named as 'the product wave function renor- 
malization group (PWFRG) method' since the RG trans- 
formation is applied to the matrix product wave func- 
tion. ^'^'^'' Numerical efficiency of the PWFRG method is 
achieved by estimating a trial wave function for the iter- 
ative calculation in the infinite system DMRG method, 
and has been confirmed through applications to classi- 



cal systems^^'^^ as well as quantum spin systems. ^^ ^^ It 
should be noted that the wave function estimation in the 
PWFRG method is of use for the finite system DMRG 
method, ^^ when the system size increase is necessary for 
preparing numerical data for several system sizes. 

In this article we report an extension of the PWFRG 
method, which can be applied to quantum systems that 
have 2 site modulation. According to this modulation, 
the RG transformation to the wave function is shifted 
by 2 lattice sites. This modified PWFRG method pro- 
vides good wave function estimation when the infinite 
system DMRG method is nearly converges to the thermo- 
dynamic limit. We also discuss how to apply the PWFRG 
method to matrix product wave functions obtained by 
the finite system DMRG method. 

In the next section we explain the matrix product 
structure of the ground state wave function. In §3 we esti- 
mate the wave function, applying the RG transformation 
to the ground state wave function. We check the numeri- 
cal efficiency of the estimated wave function in §4, where 
fidelity error in the estimation process is observed. When 
there is finite excitation gap the wave function estima- 
tion works efficiently. We discuss the estimation scheme 
proposed by McCUoch quite recently, which provides bet- 
ter estimation than the PWFRG method especially when 
the system is gapless."^^ Conclusions are summarized in 
the last section. 

2. Matrix Product Formulation 

Consider the eigenvalue problem for the ground state 
of a ID quantum system that has modulation of period 
2. An example of such systems is the dimerized 5=1/2 
Heisenberg spin chain of length 2A^, which is defined by 
the Hamiltonian 



2N-1 



h^-^) = jY: {i + 5(-i)^}s,.s,^, 



(2.1) 



where J > represents the antiferromagnetic interaction 
and where 5 the dimerization. Since we treat the MPS 
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constructed by the infinite system DMRG method, the 
state which can be further improved by the finite system 
DMRG method, we assume that the system size is even. 
The bond strength at the center, between S^ and S^ , j^, 
is J(l + 5) when TV is even and is J{1 — S) otherwise. The 
system has 2-site period even when (5 = 0, in the sense 
that total spin of the first M site J2i=i ^i alternates 
between integer when M is even and half-integer when 
odd. 

We express the ground state wave function or its vari- 
ational estimate by the notation 



■ '^N ^N ^N-1 ■ ■ 



:cri), 



(2.2) 



where a^ = ±1 (i < N) represents 2Sf in the left half of 
the system, and where a^ {i < N) represents 25'|^_|_j^_j 
in the right half. We have thus divided the whole sys- 
tem into the left and the right parts, according to the 
convention in the DMRG method. Though the system 
described by iJ^^A^) j^g^g left-right symmetry, we do not 
explicitly use it in the following formulations, in order 
not to loose generality. For example, the MPS obtained 
by the finite system DMRG method is not symmetric in 
this sense. 

Let us start from the smallest case where 2N ~ 4.^^ 
It is easy to numerically (or even manually) diago- 
nalize H^^' to obtain the ground state wave function 
\l'('*)((j]^(T2 ^2 o'l)- Since we are dealing with open bound- 
ary systems, the eigenfunctions are always real. The den- 
sity matrices for the both sides of the system 



p^ (cricT^ I CTl (72 ) 



J2 ^^^\'yy^^^^^)^^''H<^i^2^2^i) 



\a[a'^\a^a^) = Y. ^^^\<yi^2^Wi)'^^^\<^i^2^2^i) 



(2.3) 

are therefore real symmetric. Diagonalizations of p and 
p^ create block spin transformations 

P"{<^W2\^l<^2) = I]A(^2)^2W<^2l?2)^(<^l<^2le2) 

p''{a[a'^\a,a^) = y^K^2)B2{^'i^2\h)B^{^,a^\£,^) , 



I2 



(2.4) 



where A2{aia2\^2) ^^^d -^2(^i^2l'?2) ^'^^ orthogonal ma- 
trices, respectively, which represent block spin transfor- 
mations (Jia2 — > ^2 ^'^d '''i'^2 ~^ ?2- Thus the block spins 
^2 and ^2 are 4-state variables. Applying the obtained 
(faithful) block spin transformations to ^!{(j^(J2 O2 &i) we 
obtain the 'center matrix'^^ 

^2(^21^2)= 51 AK'^2l^2)^2(5-lCT2|^2)*K'^2^2^l)- 

(2.5) 

Note that the matrix A2(^2l?2) is not always diagonal, es- 
pecially when we perform the diagonalizations of density 
matrices in Eq. (2.4) independently under the condition 
that there is degeneracy in density matrix eigenvalues. 



It is possible to make A2 diagonal by applying singular 
value decomposition (SVD) directly to ^^'*', but we do 
not assume the diagonal property of the center matri- 
ces in the following. Using the obtained matrices, we can 
write ^(^' in the form of matrix product 



*('^)(CTlCr2Cr2'7l) 



(2.6) 



= ^^('^l<^2le2)A2(C2le2)S2(ai (721^2) • 
C2C2 

For a while we follow the MPS construction by the 
infinite system DMRG method. Then the next step is 
the case 2N — 6. Applying the previously obtained 
block spin transformations A2 and B2 to H^^> we ob- 
tain the super-block Hamiltonian _ff(^' that acts to the 
Hilbert space spanned by ^2? "'37 '''31 and ^2- Diagonal- 
izing H^^' we obtain the ground state wave function 
^ (^2 "^z ^z C2) i^ the renormalized linear space. In the 
same manner as we have done for Eqs. (2.3)-(2.6), we 
obtain the matrix product expression 



^^'\i2^Z-^zl2) 



(2.7) 



CaCa 



^3fe^3le3)A3(e3l4)53fe^3lC3) 



where ^3 (^2 '^31^3) and B^{^2^s,\^z) represent block spin 
transformations ^2^3 ~^ ?3 ^rid ^2 '''3 ~^ $3- The dimen- 
sion of the new center matrix A3(^3|^3) is 8. 




Fig. 1 . Matrix product representations of ground state wave func- 
tions for 2N = 4, 6, and 8. 



For convenience in the matrix product representation, 
let us introduce 1-state dummy variable ^q and ^q at the 
both ends of the system. We put them at the both ends 
of the system. For example, ^^^^ in the l.h.s. of Eq. (2.6) 
can be written as ^^''■'('?o'''i'''2 ^2 ^1 ^0); i^ ^^ neglect the 
dummy variable the original form of ^'^•' in Eq. (2.6) is 
recovered. Then how does the r.h.s of Eq. (2.6) look like? 
In order to answer this question we also introduce two 
state block spin variables ^^ and S,i, respectively, which 
is always the same as a^ and a^. Using these variables 
we define the boundary orthogonal matrices 






(2.8) 
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where d{a\b) represents Kronecker's delta (5^^. With the 
help of these boundary orthogonal matrices, we can ex- 
press ^^^' and ^(^' in the matrix product form 



vi/(4) = a.a^a^bIbI 
vi/(6) = a,a^a,a,bIbIbI 



(2.9) 



where we have changed the variables of A2 and -63 as 
^2(^10-21^2) ^'^d -62(^1^21^2); respectively. In equation 
(2.9) we regard block spin variables ^^ and ^^ as the 
matrix index and take their configuration sum, leaving 
the raw spin variables a^ and CTj . It might be better to 
regard A^ and B^ as 3-leg tensors, and r.h.s. of the above 
equation as tensor products.'^'* 

It is straight forward to extend the matrix product 
expression of the ground state wave function to arbitrary 
system size 

¥^^-> = ^1^12 • • • Aj,Aj,b], . . . bIb\ (2.10) 

= AiA2...^^_i§(2^)Bl,_i...i?2^i?I, 

where configuration sum is taken for all the block spin 
variables, and where ^1/(2^) = Aj^Kj^B]^. Figure 1 shows 
the graphical representation of vI/(2^) for 2N = 4, 6, and 
8, where cross marks represent the dummy variables ^q 
and ^0 &t the both ends, black squares the block spin 
variables, and circles the raw spin variables. From the 
computational view point, it is impossible to keep all 
the degrees of freedom in block spin transformation for 
arbitrary large system size, therefore the number of state 
of the block spin variables ^^ and f ^ are restricted at most 
m states. When there is a cut off in this sense, the r.h.s. 
of Eq. (2.10) is a variational approximation for the l.h.s. 
For example, Aj^K^B]^ is an approximation of vl/(2^) 
when the matrix dimension of A^ is restricted. 

Though we have created the matrix product wave func- 
tion in Eq. (2.10) by way of the infinite system DMRG 
method, we do not restrict ourselves about the way of 
creation of MPS in the following formulation. For ex- 
ample, we also deal MPS obtained by the finite system 
DMRG method, where the sweeping is stopped at the 
center of the system. Strictly speaking, the matrices A^ 
and Bj determined by the finite system DMRG method 
is dependent to the system size 2iV, there fore we have 
to put the system size to the matrix labels as Ai ' and 
Bl for distinction. But the notation is rather com- 
plicated, and therefore we drop the label {2N) in the 
following equations. 

Let us observe the renormalized wave function ^y^^\ 
which corresponds to the lowest energy state of the 
super-block Hamiltonian H^'^^\ It is possible to obtain 
^(^^) applying block spin transformations A^ . . . ^^v-i 
and Bl . . . Bj^_i successively to v]>(2^) as 

(2.11) 
where we have identified the wave function vpf^^) as 
a 3-leg tensor, which has (dummy) matrix indices ^q 
and ^Q in addition to the row spin variables {a} — 
cri...crjvCTjv.a-i. 



3. Wave Function Renormalization 

Suppose we have matrix product expressions for ^'^^ 
and ^e^) in Eq. (2.9), and need to obtain that of ^^(s). 
This need is fulfilled if we diagonalize the Hamiltonian 
H^^> via eigen solver such as the Lanczos method. Un- 
der the situation it is important to prepare a good trial 
(or initial) wave function for the numerical diagonaliza- 
tion process. An answer to this problem can be obtained 
from observation on the bare Hamiltonians H^'^' and 
H^^\ Since these two Hamiltonians has the same bond 
strength at the center of the system, \I'(^^ can be used as 
a trial (or variational) wave function for H^^' if we put 
two additional spins to the both ends. This construction 
is represented as 

(criCr2<^3'^4 ^4 0-302^1) = ^^^^^ ((73^4 CT4 CTg) (3.1) 



* 



(8) 
trial 



apart from the normalization factor, where the trial wave 

(8) — — 

function ^jj.^; is not dependent to a^, a2, (J2, and a-^. 
Such a construction of trial wave function can be gen- 

(2N-\-2) 

eralized to arbitrary system size 2iV, where ^^^iai i^ 
obtained from 1^(2^-2)^ Since this is a rough estimation, 
one has to improve the trial wave function afterward. 




Fig. 2. Graphical representation of the trial wave function esti- 
mation for 2A'' = 8. 



Consider the way of expressing the wave function es- 
timation in Eq. (3.1) in the renormalized subspacc. Ap- 
plying block spin transformations, which are already ob- 
tained up to 2A^ = 6, to the estimated wave function 
'^U^i{o'iO'20'3cr^ 0-4 (T3 (T2 (Ti), we obtain the renormalized 
form of the trial wave function 



^ilMs'^i^^Q 






(3.2) 

as we have done in Eq. (2.11). Figure 2 shows the graph- 
ical representation of the above wave function renormal- 
ization process applied to '^\^l^i, where we draw ^'^^' by 
its matrix product expression. In order to write Eq. (3.2) 
more transparently, we introduce dummy matrices 

^(Cl^2lCo) = I/V2 

^-i(C-2'^iie-i) = 1/V2 

BoiLi^2\Q = I/V2 
B-i{L2^i\Li) - I/V2, (3.3) 

where ^_27 C-ii Cc ^-21 C-^ij and ^q are 1-state dummy 

,(8) 
" trial 



variables. Then the ^tj-^i in Eq- (3-1) can be written as 



* 



(8) 



?t nt nt nt 



trial = A_iA,AiA2A2BlBlBlBli 

= A_iA,Ai¥''^bIbIbIi, (3.4) 
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where ^'^^H^icr^a^^i) is nothing but "ii^^^cria^a^ai) where L4 and i?4 are defined as follows 

since ^i = a^ and ^^ = ai by the definition of Ai and Bi 

in Eq. (2.8). It should be noted that in Eqs. (3.3) and 

(3.4) the matrix labels are shifted by 2, in the sense that 

Ai, A2, Bi, and B2, respectively, contains (T3, a^, Sg, and 

(J4. Substituting Eq. (3.4) into Eq. (3.2) we obtain 



^4(^41^2) - E^i^3^2 



* 



(8) 
trial 



^ aIaIa\a_,a,a,¥''^bIbIbUb,B2B, 



= E M^3'^4\^4)Lsi^3\^,)A2{^1^4\^2) 



CTj^lT^CTg CTgCTjO-j^ 



R^iLlQ = E ^1^352 



= E ^sfel^l) *^'Hei'^4 ^4 Cl) i?3(4lCl) 

where the matrix ^3(^3!^!) is defined as 

-^3(^3^1) — 2^ ^3^2^1^-1^0^1 



(3.5) 



= E B^i^3^^\^4)R3i^3\^l)B2i^l^4\Q ■ 

(3.9) 
These recursive relations in Ljy and Rj^ in the above 



(3.6) equation was first obtained empirically and has been used 

for the numerical study by use of the PWFRG method 

when it is applied to S* = 1/2 quantum spin chains. ^''^^^ 

= V" ^2('''i'''2lC2)^3(^2'''3lC3)^i(^o'''3l4) Figure 4 show the graphical representation of the relation 

cTj^a-^o-^^^ between Lg and L4, and also between Jig and i?4. The 

process of wave function estimation is drawn in Fig. 5. 



and ^3(^31^1) is in the same manner 

^3felCi) = E bIbIbIb_,b,b, 



(3.7) 



= E B2i^i^2\Q B,{^2^,\QB,{^,a,\^,) 

Figure 3 shows the graphical representation of L and 
i?, the matrices which have a function of 'adjusting' the 
dimension of block spin variables. 






Fig. 4. Construction of L^ (upper) and R^ (lower) in Eq. (3.9). 



Fig. 3. Graphical representations of L (left) and R (right) in 
Eq. (3.6). 



As ^t^jai '^^^ ^^ used as a trial wave function for 
i?*-*-*, ^triai thus obtained would be a trial wave func- 
tion for the super-block Hamiltonian H^ "> . Now we can 
perform Lanczos diagonalization of B^^' rapidly starting 
from ^tr/ai- Let us assume that we obtain the improved 
*(*)(^3cr4a-4^3) in this way. From the calculated ¥^^ 
we obtain A(C3cr4lC4)> B ^{1^^^ ^\l^ , and A4(^4|^4) as we 
have done in previous steps. 

We can extend the way of initial wave function esti- 
mation to the case 2iV = 10, where ^triai(C4'''5 ''■5C4) is 
required. This prediction is performed as 



^ trial - ^4^ ^4 



/(6)Rt 



(3.8) 



E 

?2?2 



i^4{i4\i2)^^''\i2^^-^d2)R4{Ul2 



Q- 



O 



new 

D D 



^^23^ 



Fig. 5. Graphical representation of the wave function estimation. 



It is straight forward to extend the relation in 
Eqs. (3.8) and (3.9) for arbitrary system size. This is the 
way of wave function estimation, which we call as the 
2-site shift PWFRG method. We can obtain *[riai^^' 
if we have matrix product expressions of both ^^"^^^ 
and vI'(2Af-2)^ rpj^g wave function estimation in Eq. (3.8) 
is performed using \I'(2A'-2) dji-gctly^ instead of its ma- 
trix product decomposition Aj^_^Aj^_^Bj^_-^ where ba- 
sis state restriction is imposed on A^_j^. Thus the es- 
timation becomes exact in the thermodynamic limit, 
where the matrix product wave function is position in- 
dependent. In this sense the way of estimation explained 
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here is better than the estimation in the previous for- 
mulation of the PWFRG method, ^■^'^'^ which uses trun- 
cated Aj^_^Aj^_^Blf_^, when the infinite system DMRG 
method is nearly converged. ^^ 

It should be noted that there is no need that \I>(2^) and 
^^{2N~2) ]-^g^yg ^j^g same matrices in common. For exam- 
ple, the estimation for ^t,.jjj] can be performed, if we 
have optimized ground states for {2N — 2)- and 2A^-site 
systems independently by use of the finite system DMRG 
method. In such a case the matrices A^ and B^ becomes 
system size dependent, as we have seen at the end of the 
last section. The definition of L, and R^ should be mod- 
ified according to the dependence, where the extension 
is straight forward. 

4. Convergence to the Thermodynamic Limit 

The estimated wave function 

l^Af+l---^l) 
^¥^''-^H<^3---<^N-l^N-l---^l) (4.1) 

is normally not accurate enough, since the estimated 
wave function is independent of 2 spins at each end of the 
extended system. Therefore the estimated renormalized 
wave function 

(4.2) 



,{2N+2), 
' trial 






^iZ-"'' = L,¥^''-''R^, 



might not be a good starting point for the Lanczos diag- 
onalization of i/(2^+2) . Let us check the efficiency in the 
estimation quantitatively by use of the fidelity error'^^ 



1- 



E 



T,(2W+2) 



* 



trial 



^N^ 



N"-" N + 1 '-'N + 



iCiv 



* 



(2N+2) 



i^N'^N+l '^N+1 ^n) 



i^N'^N+l f^N+l Cat) 7 (4.3) 



between normalized ^[^-^i^i and ^(2Af+2)_ -y^g observe 
the error when the wave function estimation is imple- 
mented in the infinite system DMRG method. The com- 
putational algorithm we have used for this check is as 
follows. 

(a) Diagonafizc iJ^'*) and obtain ^1/(4) = ^(4)^ ^^^ ^nd 

B2- 

(b) Create H^^^ by applying A2 and B2 to iJ^^) . Diago- 
nafizc H^^^ and obtain ^e^), A3, and B^. 

(c) Contracting A^ and B^ as Eqs. (3.6) and (3.7), re- 
spectively, to obtain L3 and R^. Set iV = 3. 

(d) Obtain ^[,^i^+^^ by applying L^ and i?]^ to §(2^-2) 
as shown in Eqs. (3.5) and (3.8). 

(e) Create the super-block Hamiltonian ij(2W+2)^ 

(f) Obtain minimum eigenvalue of i/(2A'+2) g^-^^j corre- 
sponding wave function, starting from ^^^ja^i • 

(g) Obtain Ajv-t-i ^^'^ ^N+i- By use of these transfor- 



mations, create L 



N+l 



E. 



A 



1 J^ atA^ 



N+l-^N^N-1 



and 



R 



N+l 



E^.,,,5Ui^ivSjv-iasEq. (3.9) 



(h) Set TV -Fl ^ TV and go to the step (d). 
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Fig. 6. The fidelity error in Eq. (4.3) calculated for the uniform 
S = 1/2 Heisenberg spin chain 5 = and the dimerized one 
(5 = 0.1. 



(5 = and (5 = 0.1 under the condition m = 64. In both 
cases the error decreases with the system size, and the 
decay is more rapid when S = 0.1 than 6 = 0. The be- 
havior can be explained by the fact that the effect of 
system boundary becomes weak in large size systems, 
and the center of the system is 'effectively decoupled' 
from the system boundary if the system size exceeds sev- 
eral time larger than the correlation length. As we have 
stated in the last section, the estimation becomes exact 



(2JV+2) 
trial 



^,{2N+2) jj^ ^]-jg thermodynamic limit N -^ 00, 
where the fidelity error becomes zero. 

It has been known that the DMRG method applied 
to gapless systems introduces an artificial correlation 
length, as an effect of basis state restriction to m. There- 
fore the convergence of the fidelity error with respect to 
the system size 2N is slow exponential when the system 
size 2iV is sufficiently large. In such a case, it is better to 
increase the system size as fast as possible. The recursion 
relation 



^N+l 



R 



N+l 



E 
E 



A 



I ^ M-^l' 



N+l^N-^N-1 



-Sjv+i-Rat-Stv- 



(4.4) 



can be regarded as linear transformations to Ljy and 
i?^, which have their fixed points in the limit A^ — > 00. If 
the number of block spin states does not change during 
this extension process, one can estimate these fixed 
point easily. But the number of states of the block spin 
variables are not always the same. A way of overcoming 
this difficulty is to modify the diagonalization step (f) 
as follows.^'*' ^^ 



(f) Improve the estimated wave function '^\^^g^ by 
applying the Lanczos step only once. Use the 'im- 
proved' wave function ^fimprovod ^^^ ^^^ succeeding 
processes. 



Figure 6 shows the fidelity error of the 5=1/2 Heisen- 
berg spin chain with respect to the system size 2A^ when 



When the system is gapless, the efficiency of the 
PWFRG method decreases. Quite recently McClloch re- 
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ported a new estimation scheme, which works even for 
gapless systems. Let us observe his method from the view 
angle of the wave function renormahzation. The starting 
point is to interpret the wave function as a matrix 



v&(2^)(ai...a^|S^...aO. 



(4.5) 



This 'wave function matrix' satisfies the identity relation 

where (^[/(^a'))" jg matrix inverse of v]/(2A'). McCUoch's 
way of wave function estimation is obtained by decreas- 
ing the system size of this inverse matrix by 2 

$(-)(^(--.))-$(^-) = ^(-+2), (4.7) 

where we $l ^^^ *^r ^^^ rectangular matrices 

obtained by shifting the left-right division of the system 
by 1 site. This construction is similar to the extension of 
corner transfer matrix, which has been applied to two- 
dimensional classical lattice models.^ 

It is easy to see that McCUoch's way of wave func- 
tion estimation can be performed by use of ^(^^-2) ^^^ 
vj/CsAf) iiiQ^i are created independently by the finite sys- 
tem DMRG method. Let us express ^(2^-2) and ^(^w) 
as 



vl/(2^) = A\...A'j,A'j,B'l 



N^ N 



.B'l 



(4.9) 



where Aj = A] and B^ ~ B\ are not always 



the same as A' , 



A 



(2N) 



and B', 



d(2JV) 

Bl , respec- 



tively. Substituting these matrix product wave functions 
to Eq. (4.7) we obtain 

,-1 



* 



(2Af+2) _ 



trial 



= A\...A\K\B%R^A^_,)-'L 



A'atA^S'Iv • • • S'l , (4.10) 

where the matrices L and R are defined as follows 
L = 2_^ Aj^_-^^. . .A^A I. . .A ^_i 



R 



E 

•••o-jv-i 



B 



N-1 ■ ■ 



B\B\...B' 



N-l ■ 



(4.11) 



Note that the matrices L and R becomes identity ones 
when both v]/(2^-2) ,^^^ \^(2N) ^^.^ ^j-gated succeedingly 

by the infinite system DMRG algorithm. 

5. Conclusions 

We have formulated a way of applying the PWFRG 
method for quantum spin systems which have 2-site mod- 
ulation. In order to estimate the initial wave function, we 
shift the application of renormahzation group transfor- 
mation to the wave function by 2 lattice cites. As a result, 
we obtain a recursive relation among renormalized wave 



functions. Numerical efficiency of the wave function esti- 
mation is confirmed when the method is applied to the 
dimerized S = 1/2 Heisenberg spin chain. We give an 
interpretation to McGUoch's way of wave function esti- 
mation, from the view point of wave function rcnormal- 
ization. 
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and Culture. 
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